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lifndef GGCONSTANTSJ 
♦define GGCONSTANTSJ 



const double ggFourPi ■- 12.566371; /♦ nee<)s more digits '/ 

canst double ggTwcPi = 6. 283185 J; /* needs more digits */ 

const double ggPi = 3.U15926535B919J23846j 

const double ggHalfPi = 1.57079532619489661923; 

const double ggThirdPi = 1.0471976; /* needs iwre digits V 

const double ggQuarterPi = 0.78539816; /* needs more digits V 

const double gglnversePi : 0.31930989; 

const double ggSqctl^<o = l.U4313S62373D950<88; 

const double gglnverseSqrtTvo = 0.7C7 10678; 

const double ggSqrCThree « 1.732050G07S6B877293S; 

const double ggSqrtfive = 2,23<067977499789696<; 

const double ggE t 2.718281828459045235360287; 

const double ggRad = 57.29577951308232; 

mfdef sun 

const double gglnflnicy = l.OelO; 
lindude <float.h> 
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const double gglnfinity = DBLJAX; 
*endif 

lifndef H_PI 
fdefine H.PI ggPi 
lendif 

const double ggBigEpsilon ^ O.COOl; 
const double ggSpsilon = 0.000001; 
const double ggSnallEpsilon ^ O.OOOOOOOOl; 
const double ggTinyEpsilon = O.OQOOOOOOOOOl; 



const double ggColorRatio ^ 0.0039215686274509303; 
ler.dif 

Idefine gsMinlx.y) (x < yl ? x : y 
lidefine gcHs^lx.yl Ix > yl ; x ; y 
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/* Includes required •/ 
iitdef WINCONS 
linclude <vindov£.h> 
tendif 

{Include <GL\g].h> 
linclude <GL\glut.h> 
linclude <8tdio.h> 
linclude ciiiath.h> 
lifdef USEJISTPBM 
linclude <:ppn.h> 
lendif 

/" 

* something because o[ windows. 

Jd_epri».f(i, APPENDIX page: 

/" 

• our data structure of choice 
V 

typecef struct obj I 
/• other parameters */ 
float (Mtrlx|16|; 

/■ vlev angle */ 
float viewangle; 

/• aspect ratio "I 
float aspect; 

/• z of the camera •/ 
float it] 

/• ry of the camera •/ 
float ry; 
) Obj: 

/' hold the display lists for textures V 
typedef struct texture ( 

int texl; 

int tex2; 
) Texture; 

/" 

' our global variables 
V 

/' camera settings ♦/ 
Ob] scene; 

/• texture stuff V 
Text are def; 

Texture* current jexture = idef; 

/• track the next display list nunber '/ 
int nextDLnun « 2; 

/• stuff for lighting •/ 

float lightPosK) = (2.0, 4.0, 2.0, 0); 

flodt lightDir[<l (0, 0, 1.0, 1,01; 

float lightAwb(4) = (C,3, 0,3, 0,3, 1.0); 

float lightDiffMI = {0.6, 0.6, 0.6, l.Oli 

float lightSpecUl = (0.6, 0.6, 0,6, 1.0); 

float clipDistance = 2.14; 



int left, right, parent; 
int width, height; 
lifdef USE_NETPBN 
pixel" ppirPixelB = flj 
lendif 

Guubyte* sgiCixels s 0; 

FILE' comniandsi 

int doTalceSnapshot : 0; 

fdefine HEMISPHERE! I 

void createHeniisphere (int listNuiii, Int nuiaPts, Intgeoml; 

void drawJeftO; 

void draw.rightO; 

void Keyfunsigned char, int, int); 

yi. 

' read the frame buffer and write out a ppin file 
V 

void talteSnapshotO ( 
lifdef USEJETTPBM 

static int shotta e O; 

FILE* tile; 

char name (50); 

Int index, i, j; 

/* draw everything again */ 

draw_right(); 

glFlushll; 

drawjeftd; 

glFlushll; 

/• read the pixels from the frame buffer */ 

glReadPixel8(0, 0, width, height, GLJGB, GLJNSIOIED.byte, sglPixels) 

/♦ convert them to the ppu */ 
index : 0; 

for (i '. height - 1; i >= 0; i--) ( 
for (j = 0; ] < width; jf*) ( 
PPIUSSlGHIpptiPixelslillj], 
EgiPixelslindexj, 
egi Pixels I index -f 1|, 
sgiPixelsl Index * 2|); 

index t= 3; 



/• open a file '/ 

sprintfjname, "M.ppn', shotMuia); 

Sh0tNuin4(; 

file = fopen(nanie, 'v')\ 
I* write the ppn file */ 

ppiLwriteppjn(file, ppniPixels, width, height, 255, CI; 

/* close the file '/ 

Edoselfile); 
lendif 
) 

' Read i.T the ppn files and create display listB for a texture 

* returns the dimension of the inage 

•/ 

void readTexture (Texture* t, char* filel, char* file2) ( 
lifdef USOETPBM 
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FILE 'fpl, *fp2; 

int cols, rows, i, j, index; 
pixel ♦•rapl, **inap2; 
GLubyte nexl, *tex2; 
pixval ntaxval; 
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/• open the files •/ 
fpl = fopenlfilel. 'r-l; 
fp2 ■- Eopenlfile2, T'); 
if (Ifpll ( 

fprint£(st(3err, 'Couldn't open <8\ii', filel); 

) 

if (!Cp2| I 

f print fistderr, "Couldn't open «s\n', fileZI; 

) 

I* read the ppm files */ 

mapl s ppn.readppmlfpl. icols, irows, Sinaxvall; 

fprincfistderr. 'Is: rows td \t cols >= W\n", file), rows, cols, niaxval)i 

niap2 : ppii_readppin|Ep2, kola, irows< &inaxval); 

fprinttistderr, 'ts; rows = M \t cols = »d\n', file2, rows, cols, tiaxvalli 

/• convert then V 

texl = iialloc|8izeof(GLuby:el ' rows ' cols • 3); 
tex2 = nalloclsizeoflGLubi'tel • rows • coU ' 3); 
index = 0; 

for (i : Ci i < rows; i+t| I 
for Ij = 3; j < cols; jn) I 
/• R */ 

texllindexl = PPM.GETR(mapl|i)|jl); 
tex2|indexl = FPM_G£TRIiiap2(ilIjl); 
index Hi 

/• G •/ 

texl I index) = PPJt6BTG(xapl|il[j|)) 

tex2|index! = FPH,GETC(map21il [jll; 
index tt; 

/• e •/ 

texllindexl n PPtLGETBlmaplUllj)); 
tex2|index) = PPM.GETB(map2(llijl); 
index 



/* create the textures in the left'/ 
glutSetMindowlleft); 
/• new display list'/ 
gISewListlnextDLnuir,, GLCoHPtLE); 
t->texl = nextDLnuir,; 
neKCDLnumn; 

glTexI(nage2G(CLJEXTURE_2D, 0, 3, COlS, rows, 0, GLWB, GL_UNSIGll2D_BYTE, 
texll; 

glEndListll; 

/• new display list*/ 
glMewListlnextDlnura. GV.COHPILE); 
t->tex2 r nextDLnuitij 
nextDLnuratt; 

glTexIinage2D(GL.TEXTURE_2D, 0, J, COlS, rows, 0, GL_RCB, GL.UNSIGNED.BYTE, 
tex21; 

glEndListd; 

/• create the textures in the right*/ 
glutSetWindowt right); 



/• new display list*/ 
9lNewLi6t(t^>texl, GL.COMPILEI; 

glTexllMge2D(G[._TEXTURE_2D, 0, 3, cols, rows, 0, GL.RGB, GUNSICiNED.BnB, 
texl); 

glEtidUscO: 

/• new display llstV 
glNewLi6Clt->tex2. G1.„C0MPILE) ; 

glTexIinage2D(GLTEXTURE_2D. 0, 3, cols, rows, 0, GL_RGB. GLUNSIGNEDJYTE, 
tex2); 

glEndListl); 

free(texl); 
tre«(tex2)i 
freelitapll; 
freeliiiap2l; 
llendif 



' this will initialize the display lists for the objects 

void initiallze.objectslint argc, char«*argv) ( 
float tmpKI; 

/* read in the texture */ 
readTexture(&def, argvH), argv|2|l; 

/• create henisphere left */ 
glutSetWindowdeftl; 

createHemlsphereU, 50, GL.THimE.STRIP) ; 

/♦ create hemisphere right ♦/ 
glutSetWindowlright); 

createHeroi sphered, 50, GL_TRIANGLEJTRIP); 

/* scene */ 

scene. viewangle = 50; 

scene. tz = 0; 



void diBplayj)arent(I ( 
/• clear the screen V 

glClear(a_COLORjUFFER_Bir I a.DEPTH.aUFFER JIT) j 

) 

void draw_left() ( 
float tinp(4); 

/* clear the screen •/ 
glutSetldindoMUeft); 

glClear(GL.COLORJUFFERJIT I CL.DEPTH_BUFrER_BIT) ; 

/' adjust for scene orientation */ 
glMatrixMode(GLPROJECTION) ; 
glLoadldentityl); 

gluPerspectivelEcene. viewangle, scene. aspect, 0.1, IC.O); 
glRotateflBcene.ry, 1, 0, 0); 
glTranslateflO, 0, scene. tz); 

I* draw our models V 

glHatrixModelGLMODELVIEW); 

glP-jshMatrlxl); 
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/* now draw the semi sphere V 
SlEnaclel6l._T£XTURE_2Dl; 
g]Color3f(.5, .5, .5); 
glCallList (:urrent_texture->teitll i 
glCallLlstlHEMISPHERE); 

g:Ro:atef|lBO, 0.0, 0.0. 1.01; 
glColorStC, 1. II; 

g:callList|currenl_texture->te52); 

glCallLlEClHEHISPHEREI; 

glPopMatrixO; 
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fprintf(£tderr, "left - W, gluErrorStrlng(glGetError(|)); 

) 

void display.! eft II 



draw.lettll; 
glutSwapEuffersll; 



) 



void d:aw_righMl ( 
float tnpUl; 
float height; 

/* clear the screen •/ 
glutSetWindowlriglitl; 

glClear(G[vCOU)R_EUFF£RJIT t CLDEPTHJUFFER JITI ; 

/• setup the caT,era V 

glMatrixHodelGkPROJECTICHI; 

glLoadldentityll; 

gluPerspectivelli, scene. aspect, cUpDlfitance, 10. 01; 
glTranBlate£(0, 0, -31; 
glRotatefllS, 1, 3, 01; 
glRotateEdS, 0, 1. 01; 
glOisablelGL.TEx™E_2Dl; 

/• draw our jijxlels */ 

glM3trixHo<Je|GLJ10DELVIEW); 

glPushMdtrixO; 

/* draw a cube for the camera '/ 
glPushHatrixl); 
glLcadidentityll; 
glRotateflieO, 1, 0, 01; 
glRotacef (-scene. ry, 1, 0, 01; 
glTranslaceflO, 0, scene. tJl; 

tJiplOI = tnpdl = tnp!2| = .8; 
trop|31 - 1; 

glMaterlalfvlGL front and Ma CL^SPECULAR, tijpl; 
glMaterialf(GLFP.ONTJNDJACK, GLSHININESS, 0.01; 
glHaterialfvlGLJROMDJACK. GL^AMB1ENT_AW).DIFFUSE, tip); 
glutSolldCubei.ll; 

/' now the light '/ 

tmplOj : tmpdl = tOp|2) - 0; 

trap|3| . 1) 

glLightfv(GL_LlGHTl, GLPOSITION, tmpl; 
glLightl(GL_LIGHTl, GLSPOT.CDlOFF, scene.viewangle / 2); 
tmplOj = tnpdl = fl; tnip|2l = li tnp[3l = li 
glUghtfv(GL_LIGHTl. GLSPOT.DlRECTION, tmpl; 

/• draw a cone for the viev £rustt>i» V 



glLoadldentityll; 
height s 1 - scene.tz; 
glRotatefl-scsne.ry, I, 0, 01; 
glRotdtetHS, 0, 0, 11; 
glTranslateflO, 0, -lli 
tmpIO) tmpll] . .8; 
tnpl2l : 0; 
tfflpl3] : 1.0; 

glHaterialfv(GL_FRONT_ANDJACK, GLSPECUWR, tnpl; 
glMaterialflGLJROIfrjNDJACR, GL_SHININESS, 0.0); 
glMatarialfv(Gk.FRONT_ANDJACK, GUWBIENTJND.DIFFUSE. trap); 
glutWireCone I tan (scene.viewangle ' 3.14 / 360.0) * height, height, 
glPopMatrixO; 

/♦ draw one half oE the sphere */ 

CmplOj = tnpd) = tnp[2) = trap|3) = 0; 

glMaterlalEv(GL_PROHTJWD_BACK, GL.SPECUIAR, tmp); 

gIMaterialf(GLWNT_AND.BACK, GLSHININESS, 0.0); 

glEndblelGU.TE){TURE_2D); 

traplOl = tnpdl - tmp(2] ; trap! J) = 1.0; 

glMaterialfvlGLFRONTJNDJACK. GMMBIENTJBD.DIFFUSE, tmpl; 

g]CallLiBt(current_texture->texll; 

glCallLiBt(HEM[SPHERE|; 

/♦ draw the other half •/ 

glRotatefdBO, 0, 0, II i 

glCa ULi sc I cur rentj exture->tex2 1 ; 

glCallLlstlHEHISPHEREIi 

glPopHatrixO; 

i fprlnlf(8tdeEr, 'right - »8\n', gluErrorStriiig(glGctError()ll; 
I 

; void dlsplay.rlghtd { 
draw_rightl); 
glutSvapBuffersd; 

I 
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Ideflne H_QUIT 1 
void Select lint value) 
( 

switch I value I I 
case M_QI]IT! 

exit (01; 

break; 

I 

glutPo&tRedisplaylh 

I 

void createjienud ( 
glutCreateHenulSek'ctI; 
glutAddMenuEntryrOuit", H_QUni; 
glutAitachMenu (GLUT.RIGHT^BUTTOM ; 

I 



/• Initial Ues hading iiicdel V 
void init_leftd 



glEnablelGLDEPTH.TESTI; 
glShadeHodeKGLJHOOTHI; 



/' texture Btuff ♦/ 
gl Enable (GL_TEXmE_2D); 
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glPlxelStoreilCLJNFACULlGNMENT, Blieof (GLubyte) I; 
glTexEiivf(GL.TEXTURE_ENV, 0I^TEI{TI;RE.ENVJ10DE, GU.DECAL); 
9lTexParaiiieterf<GUEXTUf(EjD, GLTEXTURE_NRAP_S, GL.CUMP); 
gnexParajneterfiGLJEXTUREjD, GL_TEXTURLWRAP_T, GLCLAMP); 
glTexParaiietetf(GL_TEXTtJRE_2D, GUEXTURLMAGJILTER, GUEAREST); 
glTexParaneterflGtjEXWREjD, GLtextur^HIN.pilter, gljeahest); 



I 

void Ir.iLrlghtll ( 

gIEnable(GlvDEPTH_TEST); 
glShadeModellGLJMOOTH); 
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/* texture stuff •/ 
glEnable(GL_TEy,TURE.20) ; 

glPixelStorei(GLl)NPACK_ALMENT, sUeoElGLutyteDi 
glTexEnvf(GLJEXTURE„ENV, GLTEnORE_EN\'JlODE, GLJ(0DIJLATE1 ; 
glTexParameterflGIuTEXmEjD, GLJEXTURE_WRAP_S, GL.CLAMP); 
glTexParaiIieterf{GIcTEXTURE.2D, GL,TEXTURE,WRAP_T, GL_CLAMP); 
g:TexParanietetffGL_TEXTURE_2D, GL_TEXTUREJ«G.FILTER, GUEAREST); 
glTexParameterf |GL_TEXTURE_2D, 6LTEXTURE_HINJILTER, GLJEARBST) ; 

/* for blending V 
glEnable(CL.BLEND); 

glBlendFunc(GLJRC_ALPHA. GL.ONLMIHUSJRCJI-PHAI ) 

/' lighting V 

glEnablelGLLIGHTINGI; 

glEnablelGL.LIGHTO); 

glLightfv{GL_:-ICHTQ, Gl.POSmON, lightPOBl; 

glLightfvlGLLIGHTO. GLAHBIENT, lightAnb); 

glLightfvlGLLIGHTO, CkDIFFL'SE, HghtDiff); 

gJLightfvlGLLIGHTO. Gl^-SPECULAR, lightSpecl; 

9lEf.able(GLLIGHTll; 
lightAff,b|21 = 0; 

lighttohlOl ^ lightAmb HI « .Bj 
lightAir,b|3) -- 1; 

glLiohtfvlGLLIGIlTl, CLJMBIENT, UghtArab); 
glLighttvlOLLIGHTl. CLJIPFUSE, UghtDitf); 
glUghtfvlGL.LIGHTl, GLSPECULAR, liglttSpecl; 



' Called when the window ie first opened and whenever 
' the window is reconfigured (moved or resized). 
»/ 

void rayReshapedri w, int hi 
I 

/* set width and height */ 
width = w; 
height = h; 

/* define the viewport */ 

glViewport (0, 0, w, hi; 

scene. aspect = 1.0MGLfloatlw/2/[GLtloat|h) 

/• reshape the subwindows '/ 

£ir$t left •/ 
glutSetMindowlleft); 
glutReshapeWindow(w/2, h); 

/• now right V 
gluLSetHindowlright); 
glu:ReshapeWindow|w/2, hi; 
glutPositionMindow|w/2« 0); 



/* a] locate memory for the snapshot thingy ♦/ 
!f (sgiPixBlsl ( 
freelsgiPixels); 

1 

sgiPixelB s Mlloc(sizeoE(CLubyte) ♦ width • height ♦ 3) 

llfdef OSEuNBTPBH 
if IppnPlxels) { 

ppraJreearraylppmPixels, height); 

1 

ppaPixels = ppni_allocarray(wldth, height); 
lend if 
1 

/*» 

* an Idle function to do automate things 
V 

tdefine move 0 
Idefine SNAP 1 
void IdleSnapshotsO ( 
static uoveOrSnap = HOVE; 

if Imoveorsnap == HOVE) I 

I* next conaiand »/ 

Keyl'r', 0< 0); 

moveOrSnap = SNAP; 
) else I 

/' snapshot V 

KeyU3, 0, 0); 

moveOrSnap = MOVE; 

) 



* Rotate both displays 
V 

void doRotatel float amt, float x, float y, float i) ( 
glutSetWindow(left); 
glRotatcflamt, x, y, zi; 
glutSetKindow(righc); 
glRotateflamt, x, y, z|; 



/♦ 

* Keyboard handler 
*/ 

void 

Key (unsigned char key, int x, Int y) 

I 

int keyint; 

float natrixKiei; 

float Iiatrix2|l6|; 

glutSetwindowUefcif 

glMatrixMode(G:-HODELVIEW|; 

glGetPloatv(GLJlODELVIEtf_MATRIX, ir,atrixl); 

glLoadldentityl); 

glutSetWindo«i(right|; 

glKatrixModeCGLHODELVlEW); 

glCetFloatv(GU101)ELVIEM.MAT|llX, fflatrix2); 

glLoadidentityl); 

/* checlt for read coiwnand */ 
if (liey == 'r'l ( 

fBcanflcoiiiiiands, 'Id', tlteylnt); 
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key = (char I keyint; 

fprintflstderr, 'read td - ic \r\\ key, keyl; 

1 

prir,tf!"»d\n", key); 
f flush (stdout): 
swicch (key) ( 
case '!': 

/• register idle function */ 

glutSeiKindowlparent); 

glutldleFunclidleSnapshotsI; APPENDIX PAGE ( 

break; 
case 13: 

/* press enter - take snapshot*/ 

doTakeSnapshot = 1; 

break; 
case 'C: 

clipDi8tan:e -= .C2; 

break; 
case 'C: 

clipDlfitance .02; 

break; 
case 'y': 

printfCry = tf\n", scene, ryl; 

scene, ry 5; 

break; 
case 'Y': 

scene, ry +t 5; 

break; 
case '2': 

scene, tz .C2; 

break; 
case 'I': 

scene, tz .02; 

break; 
case 'a': 

scene, viewangle -» 1; 

break; 
case 'A'l 

scene. viewangle i- 1; 

break; 
case S5: 

acHotace(-5. 0,0, 0,0, 1.0); 

break; 
case 57; 

doPotatelS, 0.0, 0.0, 1,0); 

break; 
case 52: 

doRotace(-5, 0.0, 1.0, 0.0); 

break; 
case 54 1 

doRotatelS, 0.0, l.C, 0.0); 

break; 
case 56: 

doRotate(5. 1,0, 0.0, 0.0); 

break; 
case 50: 

doR0tate(-5, l.O, 0.0, 0.0); 

break; 
case '?': 

fprintflstderr, "arrows - rotate the £phere\n"); 

fprintElstderr, 'a/A vlewangle\n'|; 

fprintflstderr, "c/C adjust clip plane in right viindovVn'); 

fprintf (stderr, "i/Z camera po8ilion\n*); 

fprintElstderr, "Escape quits \n"); 



break; 

case 27; /* Esc will quit •/ 

exit(l); 

break; 
default: 

fprintflstderr, "Unbound key - td \n', key); 
break; 

) 

fprintflstderr, 'clip = U viewangle = «f zdepth : \t \n', 

clipDistance, ecene.viewangle, scene, tz); 
glutSetWindowUeft); 
glHuUMatrixflmatrixl); 
glutPostRedisplayO; 
glutSetWindowlrlght); 
glMultMacrixf(inatrix2); 
glitPostSediBplayl); 

/• check for snapshot V 
if IdoTakeSnapshot) { 

takeSnapshotO; 

doTakeSnapshot = 0; 

) 



/' 

* Main Loop 

' Open window with initial willow size, title bar, 

• RGBA display mode, and handle input events. 
V 

int naindnt argc, char** argv) 
( 

glutlnitliargc, argv); 

glutlnitDisplayMode (GLur.DOUBLE I GLUTJGba); 

/♦ create the parent window •/ 
parent = glutCreateWlndow (argv[0|l; 
width 2 512; height 256; 
glut Reshapewindow (width, height); 
glutKeyboardFunclKeyi; 
glutReshaperunc (nyReshape) ; 
glutDisplayPunc (display jjarent); 
createjieniiO; 

/♦ left window */ 

left = glutCreateSubWlndow (parent, 0, 0, width / 2, height!; 

glucKeyboardFUnclKey); 

glutDisplayFunc(di8play_left); 

createjnenull; 

lnit>ft(»; 

/* righ: window •/ 

right = glutCreateSubWindowlparent, width / 2, 0, width / 2, height); 

glutKeyboardFunc(Key); 

glutDi6playPuiic(dl6play„right); 

create.iiienu()j 

init„right(); 

/* create objects */ 
initialize^objectfllargc, argv); 

/* open file •/ 
if large --^ i) I 

fprintf (stderr, "OpeDiiig Is for corai«nd8\n', argv[3l); 

conunsnds = foDen(argv(3), 'r'); 

) 

glctMainLoopO; 
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t include <niath.h> 
•include 'ggConsiants.h' 
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* ?his function takes a point In the unit square, and maps it to 

* a point on the unit hemisphere. 

* Copyright 1994 Kenneth Chiu 

* This code may be Ereely distributed and used for any purpose, commercial 
' or non-comrcial as long as attribution is maintained. 

V 

void 

mapldouble n, double y, double *xjtt. double *y.ret, double *2„recl I 

double XX, yy, offstc, theta, phi; 

K = 2'x - 1; 
y ^ 2»y - V, 

if (y > -X) ( /♦ Above y « -x •/ 

if (y < x( ( /• Below y s x */ 

xx i x; 

if (y > OM /* Alwve x-axis V 

/• cerr«'Octant 1 

offset * 0; 
yy yj 

1 else ( /* Below and including x-axis */ 

/' cerr«'Octant 8 "i*/ 

offset = l7»ggPi)/«; 
yy * X t y; 



/' Above and including y = x •/ 
/' Right of y-axis •/ 
/' cerr«"Octant 2 



) else ( 
XX s y; 
if (X > 0) ( 



offset ! ggPi/4; 
yy = (y - x),- 

) else ( /' Left of and including y-axis */ 

cerr«*Octant 3 ";♦/ 



offset = (ggTwoPi)/4; 
yy -x; 



I 



I else \ /♦ Below and including y s -x */ 

if (y > x! I /* Above y : x */ 

XX B -x; 

if (y > 0) I /* Above x-axis •/ 

cerr«"Octant 4 ■;*/ 

offset = (3'ggPl(/l; 
yy = -X ■ y; 

) else I /* Below and including x-axls •/ 



cerr«'Octant 5 ";V 



offset s (ggFourPil/4; 
yy = -y; 

I 

) else ( 
XX = -y; 

if Ix > 01 ( /* Right of y-axis »/ 



/♦ Below and including y = x ♦/ 



cerr« "Octant 7 ';♦/ 

offset = (6*ggPil/4i 
yy = x; 

I else ( /* Left of and including y-axis V 

if (y != 0) 1 

cerr«'Octant 6 •;♦/ 

offset e (5*ggPi)/4; 
yy = X - y; 
} else { 

cerr« 'Origin ';*/ 

•x_ret = 0; 
♦y.ret = 1; 
•z.ret = Oi 
return; 



theta = acosd - xx*xxl; 

phi B offset t (ggPi/4l*(yy/xx)( 

*x_ret = Bln(thetal*cos(phl); 

*y_ret = cos (theta); 

'2.ret s sin(thetal*8in(phll) 



' This function takes a point in the unic hemisphere, and naps it to 
' a point on the unit square, 

' Copyright 1994 Keneth Chiu and Kurt Ziianerman 



* This code may be freely distributed and used for any purpose, commercial 

♦ or non-coniniercial as long as attribution is maintained. 

♦y 

void 

uraiiap (double x, double y, double i, double 'x ret, double 'y ret) 
( 

double XX, yy, offset, theta, phi; 
theta 3 acoalyl ; 

if(theta < .0000001) /' vertical center */ 

I 

*x„ret = 0.5; 
•y.ret = 0.5; 
return; 

) 

else 



DoSTTCr OSZSciiSDU 



APPENDIX PAGE { 



double cosphi.sinphi; 
cosphi : x/Bin(theta)i 

cosphi = gg.Minlcosphi, 1.01; /* hack for now V 

cosphi = ggMaxIcosphl, -l.D); 
sini^i : z/sinltheta); 

sinphi s ggMinlBinphl, l.O); /* hack for now V 

sinphi - ggMaxIslnphi, -1.0); 

if (sinphi >= 0 ) 

phi acos (cosphi I; 
else if (siaphi < 0 W C08phii< 0) 

phi -- -acos (cosphi 1) 
else 

phi : asinisinphil; 

XJi s sqrtd - yl; 

if (phi < -0 ' KJDH) 
( 

/* cerr«'5th octant ';'/ 
yy = ((phi i * xx)/(H.pi/«); 
*x.ret : -XX! 
•y_ret s -yy; 

) 

else if (phi < -M_Pl/2l 
( 

/•cerr«"6ch octant ";V 

yy = ((phi 4 (3'M_Pl/4)) * xx)/(MJI/4l; 

•y.ret = -xx; 

•x_ret n *yj6t + yy; 

I 

else if (phi < -(iLPII/41 

( 

/•cerr«'7th octant •;*/ 

yy llphi ♦ H,Pi/2l ♦ xx)/((CPl/<); 

•x.ret = yy; 

*y_rei = -xx; 



) 

else 
( 



*y_ret = xx; 



/♦cerr«Mth octant ';V 

offset = 3*H_PI/4; 

yy = (Iphi - offset) ♦ xx|y(M.PI/4) 

•x.ret = -xxj 

'y^ret e -yy - *x,ret; 



*y-ret : (•y_ret t l)/(2 ♦ ggEpsilonl; 
Xret «Cx_ret t + ggEpsilon); 



I 

else if (phi < 0) 
( 

/•cerr«' 8th octant •;•/ 

yy = ((phi 4 M,PI/4) ♦ xx)/lH.PI/41; 

'x_ret : XX; 

*y_rBt : yy - *x_ret; 

I 

else if(pM < (H_PI!/4) 

/•cerr«"lst octant ';*/ 
yy = I (phi) • ;oti/(M_PI/4); 
•x_ret 0 XX; 
*y_ret = yy; 

) 

else if (phi < N.Pl/2) 

( 

/♦cerr«"2nd octant 

yy ((phi -UI/O * ia)/(H-PI/4|; 

*y„tet : XX; 

'x.ret t *y_rec - yy; 

j 

else iftphl < 3*!LPI/4) 

( 

/'cerr« '3rd octant •;•/ 

yy = ((phi - V^m * xxl/(HjI/4l, 

'x_ret = -yy; 



iDidiiErF'tlDi 
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/ * "/ 

/• warp.h •/ 

iifndef NARP.DOTH 
Mefine WARP.KTH 

void itBpfdouble k, double y, double 'x.ret, double *y.ret, double •z_ret|; 
void unmap (double x, double y, double z, double *x_r«, double *y_ret); 
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